Explicit solution of the optimal fluctuation problem 
for an elastic string in a random medium 
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The free-energy distribution function of an elastic string in a quenched random potential, Pl(F), is 
investigated with the help of the optimal fluctuation approach. The form of the far-right tail of Pl (F) 
is found by constructing the exact solution of the nonlinear saddle-point equations describing the 
asymptotic form of the optimal fluctuation. The solution of the problem is obtained for two different 
types of boundary conditions and for an arbitrary dimension of the imbedding space 1 + d with d 
from the interval < d < 2. The results are also applicable for the description of the far-left tail of 
the height distribution function in the stochastic growth problem described by the d-dimensional 
Kardar-Parisi-Zhang equation. 
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I. INTRODUCTION 

In many physical situations, the behavior of some ex- 
tended manifold is determined by the competition be- 
tween its internal elasticity and interaction with exter- 
nal random potential, which at all reasonable time scales 
can be treated as quenched. A number of examples of 
such a kind includes domain walls in magnetic materials, 
vortices and vortex lattices in type-II superconductors, 
dislocations and phase boundaries in crystals, as well as 
some types of biological objects. It is expected that the 
presence of a quenched disorder makes at least some as- 
pects of the behavior of such systems analogous to those 
of other systems with quenched disorder, in particular, 
spin glasses. 

Following Ref. [l], in the case when internal dimension 
of an elastic manifold is equal to 1 (that is, an object 
interacting with a random potential is an elastic string) , 
the systems of such a kind are traditionally discussed 
under the generic name of a directed polymer in a ran- 
dom medium. In such a case the problem turns out to 
be formally equivalent [2j to the problem of a stochas- 
tic growth described by the Kardar-Parisi-Zhang (KPZ) 
equation 

3] 

. which in its turn can be reduced to the Burg- 
ers equation 0] with random force (see Refs. [H and for 
reviews) . 

The investigation of Pl(F), the free-energy distribu- 
tion function for a directed polymer (of a large length L) 
in a random potential, was initiated by Kardar [7|, who 
proposed an asymptotically exact method for the calcu- 
lation of the moments Z n = Z n of the distribution of 
the partition function Z in a (1 + 1) -dimensional system 
(a string confined to a plane) with a ((-correlated ran- 
dom potential and made an attempt of expressing the 
moments of Pl(F) in terms of Z n . Although soon af- 
ter that, Medina and Kardar Q (see also Refs. [H and 
0) realized that the implementation of the latter task is 
impossible, the knowledge of Z n allowed Zhang to 
find the form of the tail of Pl{F) at large negative F. 
The two attempts of generalizing the approach of Ref. 
to other dimensions were undertaken by Zhang [llT | and 



Kolomeisky [12j |. 

Quite recently, it was understood [H, [3] that the 
method of Ref. 10 allows one to study only the most dis- 
tant part of the tail (the far-left tail), where Pl{F) is not 
obliged to have the universal form Pl (F) — P* {F/F* ) /F* 
(with F* oc L w ) it is supposed to achieve in the thermo- 
dynamic limit, L — ► oo. For (1 + 1) -dimensional sys- 
tems the full form of the universal distribution function 
is known from the ingenious exact solution of the polynu- 
clear growth (PNG) model by Prahofer and Spohn [l5l ]. 
However, there is hardly any hope of generalizing this ap- 
proach to other dimensions or to other forms of random 
potential distribution. 

One more essential step in the investigation of different 
regimes in the behavior of Pl(F) in systems of different 
dimensions has been made recently [lj, [3| on the basis 
of the optimal fluctuation approach. The original version 
of this method was introduced in the 1960s for the inves- 
tigation of the deepest part of the tail of the density of 
states of quantum particles localized in a quenched ran- 
dom potential 0, [TtI EH • Its generalization to Burgers 
problem has been constructed in Refs. Qj^ and [2(| but 
for the quantities which in terms of the directed poly- 
mer problem are of no direct interest, in contrast to the 
distribution function Pl(F) studied in Refs. [H and [Til. 
Another accomplishment of Refs. [IH and UH consists in 
extending the optimal fluctuation approach to the region 
of the universal behavior of Pl(F), where the form of this 
distribution function is determined by an effective action 
with scale-dependent renormalized parameters and does 
not depend on how the system is described at microscopic 
scales. 

In the current work, the results of Refs. [l3| and [U 
describing the behavior of Pl(F) at the largest posi- 
tive fluctuations of the free energy F (where they are 
not described by the universal distribution function) are 
rederived at a much more quantitative level by explic- 
itly finding the form of the optimal fluctuation which is 
achieved in the limit of large F. This allows us not only 
to verify the conjectures used earlier for finding the scal- 
ing behavior of S(F) = — \h[Pl(F)] in the corresponding 
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regime, but also to establish the exact value of the nu- 
merical coefficient entering the expression for S(F). For 
brevity, we call the part of the right tail of Pl (F) studied 
below the far-right tail. The outlook of the paper is as 
follows. 

In Sec. [Ill we formulate the continuous model which 
is traditionally used for the quantitative description of a 
directed polymer in a random medium and remind how 
it is related to the KPZ and Burgers problems. Sec. HTT1 
briefly describes the saddle-point problem which has to 
be solved for finding the form of the most optimal fluc- 
tuation of a random potential leading to a given value 
of F, In Sec. [iVl we construct the exact solution of the 
saddle-point equations introduced in Sec. [HI] for the case 
when the displacement of a considered elastic string is 
restricted to a plane (or, in other terms, the transverse 
dimension of the system d is equal to 1). We do this 
for sufficiently large positive fluctuations of F, when the 
form of the solution becomes basically independent on 
temperature T and therefore can be found by setting T 
to zero. 

However, the solution constructed in Sec. IIVI turns out 
to be not compatible with the required boundary condi- 
tions. Sec. [V] is devoted to describing how this solution 
has to be modified to become compatible with free initial 
condition and in Sec. IVI[ the same problem is solved for 
fixed initial condition. In both cases, we find the asymp- 
totically exact (including a numerical coefficient) expres- 
sion for S{F) for the limit of large F. In Sec. ECD the 
results of the two previous sections are generalized for 
the case of an arbitrary d from the interval < d < 2, 
whereas the concluding Sec. IVIIII is devoted to summa- 
rizing the results. 



II. MODEL 

In the main part of this work, our attention is focused 
on an elastic string whose motion is confined to a plane. 
The coordinate along the average direction of the string 
is denoted t for the reasons which will become evident 
few lines below and x is string's displacement in the per- 
pendicular direction. Such a string can be described by 
the Hamiltonian 



H{x(t)} 





~dx(t)~ 




dt 



-i 2 



+ V[t,x(t)} 



(1) 



where the first term describes the elastic energy and 
the second one the interaction with a random potential 
V(t,x), with L being the total length of a string along 
axis t. Note that the form of the first term in Eq. |T|) re- 
lies on the smallness of the angle between the string and 
its preferred direction. 

The partition function of a string which starts at t = 
and ends at the point (t, x) is then given by the func- 
tional integral which has exactly the same form as the 
Euclidean functional integral describing the motion of a 



quantum-mechanical particle (whose mass is given by J) 
in a time-dependent random potential V(t, x) (with t 
playing the role of imaginary time and temperature T of 
Plank's constant K). As a consequence, the evolution of 
this partition function with the increase in t is governed 
[2j by the imaginary-time Schrodinger equation 



Tz = 



z(t, x) 



(2) 



Here and below, a dot denotes differentiation with respect 
to t and V differentiation with respect to x. 

Naturally, z(t, x) depends also on the initial con- 
dition at t = 0. In particular, fixed initial condi- 
tion, x(t = 0) = xq, corresponds to z(0,x) = S(x — xq), 
whereas free initial condition (which implies the absence 
of any restrictions d, [HI on x at t = 0) to 



z(0,x) = 1 



(3) 



Below, the solution of the problem is found for both these 
types of initial condition. 

It follows from Eq. ([2]) that the evolution of the free 
energy corresponding to z(t,x), 



f{t,x) = -Tin [z{t,x)\ 



(4) 



with the increase in t is governed [2J by the KPZ equation 
i 



2 J 



(V/) 2 - ^V 2 / = V(t,x) 



(5) 



with the inverted sign of /, where t plays the role of time 
and v = T/2J of viscosity. On the other hand, the deriva- 
tion of Eq. (O with respect to x allows one to establish 
the equivalence [2| between the directed polymer problem 
and Burgers equation 0] with random potential force 



uWu - vv 2 u = j _1 vy(t, x) , 



(6) 



where u(t,x) — V/(£, x)j J plays the role of velocity. Note 
that in terms of the KPZ problem, the free initial con- 
dition ([3]) corresponds to starting the growth from a flat 
interface, /(0, x) — const, and in terms of the Burgers 
problem, to starting the evolution from a liquid at rest, 
u(0,x) = 0. 

To simplify an analytical treatment, the statistic of a 
random potential V(t, x) is usually assumed to be Gaus- 
sian with 



V(t,x) = 0, V(t,x)V(t',x') = U(t-t',x-x'), (7) 

where an overbar denotes the average with respect to dis- 
order. Although the analysis below is focused exclusively 
on the case of purely 5-functional correlations, 

U(t -t',x- x) = U 6(t - t')S(x - x) , (8) 

the results we obtain are applicable also in situations 
when the correlations of V(t, x) are characterized by a fi- 
nite correlation radius £ because in the considered regime, 
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the characteristic size of the optimal fluctuation grows 
with the increase in L and therefore for large-enough L, 
the finiteness of £ is of no importance and an expression 
for U(t — t',x — x') can be safely replaced by the right- 
hand side of Eq. (|SJ) with 



+ 00 



dt 



+ OO 



dx U(t, x) . 



(9) 



III. OPTIMAL FLUCTUATION APPROACH 

We want to find probability of a large positive fluctu- 
ation of free energy of a string which at t — L is fixed 
at some point x = Xl, . It is clear that in the case of free 
initial condition, the result cannot depend on xl, so for 
the simplification of notation, we assume below xl = 
and analyze fluctuations of F = f(L, 0) — /(0, 0). As in 
other situations [H, [I?], EH], tne probability of a suffi- 
ciently large fluctuation of F is determined by the most 
probable fluctuation of a random potential V (t, x) among 
those leading to the given value of F. 

In its turn, the most probable fluctuation of V(t, x) 
can be found 21| by looking for the extremum of the 
Martin-Siggia-Rose action [H, [H, [H| corresponding to 
the KPZ problem, 



S{f,V} 



1 



1, 



r+co 

, dt I dxl --V 2 
U J ]_„ I 2 







+ V 



/ + ^(v/) 2 -^v 2 / 



, (10) 



both with respect to / = f(t, x) and to a random poten- 
tial realization V = V(t, x). The form of Eq. (flT)]) ensures 
that its variation with respect to V(t, x) reproduces the 
KPZ equation (0, whose substitution back into Eq. HT 
reduces it to the expression 



S{V} 



1 

2U ( 



dt 



Jo 



dxV 2 {t,x) . 



(11) 



determining the probability of a given realization of a 
random potential, "P{V} oc exp(—S{V}). On the other 
hand, variation of Eq. (fT0|) with respect to f(t,x) shows 
that the time evolution of the optimal fluctuation of a 
random potential is governed by equation [25[ 

V + —V (Wf) + vV 2 V = 0, (12) 

whose form implies that the integral of V(t, x) over dx is 
a conserved quantity. 

Our aim consist in finding the solution of Eqs. ([5]) and 
([12]) satisfying condition 



/(L,0)-/(0,0) = F, 



(13) 



as well as an appropriate initial condition at t = 0. The 
application of this procedure corresponds to calculating 



the full functional integral determining Pl(F) with the 
help of the saddle-point approximation. In the framework 
of this approximation, the condition (|13[) [which formally 
can be imposed by including into the functional integral 
determining Pl(F) the corresponding ^-functional fac- 
tor] leads to the appearance of the condition on V(t,x) 
at t = L [H, 



V(L,x) = X6(x), 



(14) 



where, however, the value of A should be chosen to ensure 
the fulfillment of condition (JT3J) . 

The conditions for the applicability of the saddle-point 
approximation for the analysis of the far-right tail of 
P L {F) are given by S > 1 and F > JU^L/T 4 . The 
origin of the former inequality is evident, whereas the 
fulfillment of the latter one ensures the possibility to ne- 
glect the renormalization of the parameters of the sys- 
tem by small-scale fluctuations [26| . We also assume that 
F ^> T, which ensures that the characteristic length scale 
of the optimal fluctuation is sufficiently large to ne glec t 
the presence of viscous terms in Eqs. (0) and (fT2j) [2q |. 
This allows us to replace Eqs. © and (fT2")) by 



/ + ^(V/) 2 = V, 
V+jV(VVf) - 0, 



(15a) 
(15b) 



which formally corresponds to considering the original 
(directed polymer) problem at zero temperature, T = 0, 
where the free energy of a string is reduced to its ground 
state energy. In accordance with that, in the T = limit 
f(L,0) is given by the minimum of Hamiltonian ((T|) on 
all string's configurations x(t) which at t = satisfy a 
chosen initial condition and at t = L end up at x(L) = 0. 

Exactly like Eq. JHJ), Eq. (|15bj) implies that V(t,x) 
behaves itself like a density of a conserved quantity, but 
takes into account only the nondissipative component to 
the flow of V given by Vu, where 



u(t,x) = J- l Vf(t,x) 



(16) 



plays the role of velocity. Naturally, for v = the time 
evolution of u is governed by the nondissipative version 
of the force-driven Burgers equation ©, 

u + uVu = VV/J . (17) 



IV. EXACT SOLUTION OF THE 
SADDLE-POINT EQUATIONS 

It is clear from symmetry that in the optimal fluctua- 
tion we are looking for, both f(t,x) and V(t,x) have to 
be even functions of x. After expanding them at x = in 
Taylor series, it is easy to verify that an exact solution of 
Eqs. (fT5]) can be constructed by keeping in each of these 
expansions only the first two terms: 



f(t,x) = J [A(t) — B(t)x 2 ] 
V{t,x) = J [C(t) - D(t)x 2 ] 



(18a) 
(18b) 
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Substitution of Eqs. (fT8|) into Eqs. (jT5|) gives then a 
closed system of four equations, 



A = C, 

B = 2B 2 + 

C = 2BC , 

D = 6BD, 



D. 



(19a) 
(19b) 
(20a) 
(20b) 



which determines the evolution of coefficients A, B, C 
and D with the increase in t. 

It is easy to see that with the help of Eq. (|20b[) , D(t) 
can be expressed in terms of B(t), which allows one to 
transform Eq. (|19b[) into a closed equation for Bit), 



6 / dt'B{t') 

to 



(21) 



B = 2B 2 + L>(t )exp 
After making a replacement 

(22) 

Eq. pT|) is reduced to an equation of the Newton's type, 

(23) 



£ = 0, 



where a = 2D(t)0 3 (t) is an integral of motion which does 
not depend on t. 

Eq. (|23]) can be easily integrated which allows one 
to ascertain that its general solution can be written as 
4>{t) = (f>o$[(t — to)/L*\, where to and 4>o = </>(to) are ar- 
bitrary constants, 

L - = wk^ m 

plays the role of the characteristic time scale, and $(rj) 
is an even function of its argument implicitly defined in 
the interval — 1 < r/ < 1 by equation 

- $) + arccos\/¥ = . (25) 

With the increase of \rj\ from to 1, $(77) monotoni- 
cally decreases from 1 to 0. In particular, on approaching 
r\ = ±1, the behavior of <&(??) is given by 



$(r,)«[(3V4)(l~H)] 2/3 



(26) 



Since it is clear from the form of Eq. (|22|) that the con- 
stant 4>o drops out from the expression for -B(t), one with- 
out the loss of generality can set 0o = 1 and 



t-t 



(27) 



The functions A(t), B(t), C(t) and D(t) can be then 
expressed in terms of (j> = (j>{t) as 

Co 



A(t) = A Q + sign(t - t ) 



1/2 



B{t) = sign(t-to) [r>o(l 
C(t) = Co/0, 
D(t) = Co/0 3 , 



3lV2 



\R> , (28a) 

(28b) 
(28c) 
(28d) 



where A = A(t ), C = C{t ) and D a = D(t ). 

Thus we have found an exact solution of Eqs. (fT5|) in 
which f(t,x) is maximal at x — (for t > to) and the 
value of /(t, 0) monotonically grows with the increase in 
t. However, the optimal fluctuation also have to satisfy 
particular boundary conditions. The modifications of the 
solution (|2~T|) - (f2"5|) compatible with two different types of 
initial conditions, free and fixed, are constructed in Sees. 
IVland lVll respectively. 



V. FREE INITIAL CONDITION 

When the initial end point of a polymer (at t — 0) 
is not fixed (that is, is free to fluctuate), the boundary 
condition at t = can be written as z(0, x) = 1 or 



/(Q,x) = 0. 



Apparently, this condition is compatible with Eq. (|18a|) 
and in terms of functions A(t) and B(t) corresponds to 



A(0)=0, B(0) = 0, 



(29) 



from where 0(0) = and to = 0. However, it is clear that 
the solution described by Eqs. (f2"T)) - (|2"9"|) cannot be the 
optimal one because it does not respect condition (p~4|) 
which has to be fulfilled at t = L. Moreover, this solution 
corresponds to an infinite action and the divergence of 
the action is coming from the regions where potential 
V(t, x) is negative, which evidently cannot be helpful for 
the creation of a large positive fluctuation of /(t,0). 

From the form of Eqs. (Q]) and (fTTj) . it is clear that any 
region where V(t, x) < cannot increase the energy of a 
string but makes a positive contribution to the action. 
Therefore, in a really optimal fluctuation with F > 0, 
potential V(t, x) should be either positive or zero. In par- 
ticular, since just the elastic energy of any configuration 
x(t) which somewhere crosses or touches the line 



,{t) = [2F{L-t)/J] 



1/2 



(30) 



and at t = L ends up at x(L) = is already larger than 
F, there is absolutely no reason for V(t, x) to be nonzero 
at least for \x\ > x*(t). 

It turns out that the exact solution of the saddle-point 
equations (|15[) in which potential V(t, x) satisfies bound- 
ary condition (fl~4|) and constraint V(t, x) > can be 
constructed on the basis of the solution found in Sec. IIVI 
just by cutting the dependences (JTHJ) at the points 



x = ±X(t), X(t) 



\c(ty 


1/2 

= 


'OA 


[D(t)_ 







1/2 

<t>{t): 



(31) 



where V(t, x) is equal to zero, and replacing them at 
I a; I > X(t) by a more trivial solution of the same equa- 
tions with V(t, x) = which at x = ±A(t) has the same 
values of f{t 7 x) and u(t,x) as the solution at \x\ < X. 
Such a replacement can be done because the flow of V 
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through the moving point x = X(t) in both solutions is 
equal to zero. In accordance with that, the integral of 
V(t, x) over the interval —X(t) < x < X(t) does not de- 
pend on t. It is clear from Eq. (j3"Tj) that X(t) is maximal 
at t = to and at t > to monotonically decreases with the 
increase of t. 

The form of fit, x) at |x| > X(t) is then given by 



f(t,x)=f[t,X(t)] 



\x\ 

J I dx' u (t,x') 
lx(t) 



(32) 



where Uo(t, x) is the solution of Eq. (fT7|) with zero right- 
hand side in the region x > X(t) which at x = X(t) 
satisfies boundary condition 



u [t,X(t)]=v(t). 



(33) 



In Eq. (|33[) . we have taken into account that in the solu- 
tion constructed in Sec. HV1 u[t. X(i)] = -2B(t)X(t) co- 
incides with 



v(t) 



dX 
~dt 



1/2 



= -2VCo(<^- 1 -l), (34) 



the velocity of the point x — X(t). This immediately fol- 
lows from Eq. (|2"2"|) and ensures that the points where 
spacial derivatives of u(t, x) and V(t,x) have jumps al- 
ways coincide with each other. It is clear from Eq. ([34]) 
that v(to) = 0, whereas at t > to, the absolute value of 
v (t) < monotonically grows with the increase in t. 

Since Eq. (|17|) with vanishing right-hand side implies 
that the velocity of any Lagrangian particle does not de- 
pend on time, its solution satisfying boundary condition 
(l33l) can be written as 



u (t,x) = v[r(t,x)} , 



(35) 



where function r(t, x) is implicitly defined by equation 



x = X(t) + (t - t)v(t) 



(36) 



Monotonic decrease of v(t) < with the increase in 
r ensures that in the interval X(t) < x < Xo = -^(^o), 
Eq. (|36p has a well-defined and unique solution which 
at fixed t monotonically decreases from t at x = X(t) 
to at x = Xo- In accordance with that, uo{t,x) as a 
function of x monotonically increases from v(t) < at 
x = X(t) to at x = Xo- For free initial condition (im- 
plying to = 0), the form of the solution at x > Xq remains 
the same as in the absence of optimal fluctuation, that 
is, uo[t,x > Xo] = 0. The fulfillment of the inequality 
duo(t,x)/dx < in the interval x > X(t) demonstrates 
the absence of any reasons for the formation of additional 
singularities (such as shocks) , which confirms the validity 
of our assumption that the form of the solution can be 
understood without taking into account viscous terms in 
saddle-point equations (J5]) and (fT2|l . 

Substitution of Eqs. and ([35]) into Eq. <JH3> and 
application of Eqs. (|15a|) and (JTTJ) allow one to reduce 



Eq. d32|) to 



r(t,x) 



/(*,*) = - / dr'(t-r') 



,,dv 2 {T') 



dr> 



(37) 



from where it is immediately clear that on approaching 
x = Xo, where r(t,x) tends to zero, f(t,x) also tends 
to zero, so that at \x\ > X , the free energy f{t,x) is 
equal to zero (that is, remains exactly the same as in the 
absence of optimal fluctuation). However, for our pur- 
poses, the exact form of the solution at |x| > X(t) is of 
no particular importance, because this region does not 
contribute anything to the action. 

It is clear that the compatibility of the constructed 
solution with condition (fl"4|) is achieved when the inter- 
val [— X{t),X{t)] where the potential is non-vanishing 
shrinks to a point. This happens when the argument of 
function $ in Eq. (|2"T|) is equal to 1, that is when 



which for to — corresponds to L„ = L and 



(38) 



(39) 



On the other hand, Eq. (|28a[) with Ao = gives 
A(L) = (tt/2)C /D 1 /2 = 2LC . With the help of the con- 
dition A(L) = F/J, following from Eq. ((T3J) this allows 
one to conclude that 



C = 



F 
2JL 



(40) 



Thus, for free initial condition, the half width of the re- 
gion where the optimal fluctuation of a random potential 
is localized is equal to 



9i 

Do 



1/2 



2 f2FL\ 



1/2 



TT \ J 



(41) 



at t = (when it is maximal) and monotonically de- 
creases to zero as X(t) = X Q(t/L) when t increases to 
L. On the other hand, V(t, 0), the amplitude of the po- 
tential, is minimal at t = (when it is equal to F/2L) 
and monotonically increases to infinity. In the beginning 
of this section, we have argued that V(t, x) has to van- 
ish at least for \x\ > x*(t) = [2F(L - ^/J] 1 / 2 and indeed 
it can be checked that X(t) < x*(t) at all t, the maxi- 
mum of the ratio X(t)/x„(t) being approximately equal 
to 0.765. 

In the case of free initial condition, the optimal fluc- 
tuation of a random potential at T = has to ensure 
that E(xo), the minimum of H{x(t)} for all string's con- 
figurations with x(0) — xq and x{L) — 0, for all values 
of xq should be equal or larger than F. In particular, 
for any xo from the interval \xq\ < X where the poten- 
tial is nonzero, the corresponding energy E(xo) has to be 
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exactly equal to F, otherwise there would exist a possi- 
bility to locally decrease the potential without violating 
the condition E(xq) > F. 

The configuration of a string, x(t), which minimizes 
H{x(t)} in the given realization of a random potential 
for the given values of xo and x{L), at < x < L has to 
satisfy equation 



J 



d 2 x dV(t,x) 



dt 2 



dx 



0. 



(42) 



which is obtained by the variation of Hamiltonian ([I]) 
with respect to x. It is not hard to check that for the opti- 
mal fluctuation found above, the solution of this equation 
for an arbitrary xq from the interval — Xq < xq < Xq can 
be written as 



x(t) 



Xn 



X(t). 



(43) 



All these solutions have the same energy, E(xo) = F. 

The value of the action corresponding to the optimal 
fluctuation can be then found by substituting Eqs. (|18bp . 
(l2"8"c)) and (|28d|) into the functional (JTTJ) , where the in- 
tegration over dx should be restricted to the interval 
— X(t) < x < X(t), which gives 



SI 



free 



8 C 5/2 J 2 

15 D y 2 Uo 



dt 

W) 



4nC 5 /2 J 2 
15 D Uo 



(44) 



The integral over dt in Eq. (|44|) can be calculated with 
the help of replacement dt/cj> — d<p/(4>4>) and is equal to 
tt/2D . Substitution of relations ([39]) and (|40[) allows 
one to rewrite Eq. (j44[) in terms of the parameters of the 
original system as 



S bcc (F,L) = K 



F 5/2 



C/ J 1/2 £ 1/2 



K = 



8V2 
157r 



(45) 



The exponents entering Eq. (|45[) have been earlier 
found in Ref. [TH from the scaling arguments based on 
the assumption that for large L, the form of the opti- 
mal fluctuation involves a single relevant characteristic 
length scale with the dimension of x which algebraically 
depends on the parameters of the system (including L) 
and grows with the increase of L [27|. The analysis of 
this section has explicitly confirmed this assumption and 
has allowed us to find the exact value of the numerical 
coefficient K. 

Since the characteristic length scale of the solution we 
constructed is given by Xo = X(t ) ~ (FL/J) 1 / 2 , the ne- 
glect of viscosity v remains justified as long as the charac- 
teristic relaxation time corresponding to this length scale 
T rc i ~ X 2 /v ~ FL/T is much larger than the time scale 
of this solution L, which corresponds to 



F> T. 



(46) 



Another condition for the validity of Eq. (j4"5"|) is the con- 
dition for the direct applicability of the optimal fluctu- 
ation approach. One can disregard any renormalization 



effects as long as the characteristic velocity inside opti- 
mal fluctuation is much larger [l4[ than the characteristic 
velocity of equilibrium thermal fluctuations at the length 
scale x c ~ T 3 /JUo, the only characteristic length scale 
with the dimension of x which exists in the problem with 
(5-functional correlations (that is, can be constructed from 
T, J and Uo). In terms of F, this condition reads 



(47) 



F > U^JL/T 4 



It is easy to check that the fulfillment of conditions (|46|) 
and (|4"T|) automatically ensures S > 1, which also is a 
necessary condition for the applicability of the saddle- 
point approximation. 

For L 3> L c , where L c ~ T 5 /JUq is the only charac- 
teristic length scale with the dimension of L which exists 
in the problem with (5-functional correlations, condition 
(|46|) automatically follows from condition (|47|) which can 
be rewritten as F > (L/L C )T. Thus, for a sufficiently 
long string (with L 3> L c ), the only relevant restriction 
on F is given by Eq. (|4"T|) . 



VI. FIXED INITIAL CONDITION 

When both end points of a string are fixed 
[x(0) = Xq, x(L) = xl], one without the loss of gener- 
ality can consider the problem with xo = Xj,. Due 
to the existence of so-called tilting symmetry [28]], the 
only difference between the problems with xo — xl and 
io/il consists in the shift of the argument of Pl (F) 
by AF = J(xl — xq) 2 /2L. For this reason, we consider 
below only the case xo — xl = 0. 

When a string is fastened at t = to the point x = 0, in 
terms of z(t,x), the boundary condition at t = can be 
written as z(0,x) oc 6(x). In such a case, the behavior of 
f(t, x) at t — > is dominated by the elastic contribution 
to energy, which allows one to formulate the boundary 
condition in terms of /(£, x) as fbl ] 



lim 

t->o 



f(t,x)-f<»(t,x) 



= 0, 



(48) 



where f^(t,x) — Jx 2 /2t is the free energy of the same 
system in the absence of a disorder. Since we are explic- 
itly analyzing only the T — ► limit, we omit the linear 
in T contribution to the expression for f (Q) (t,x) which 
vanishes in this limit. The fulfillment of condition ([48]) 
can be ensured, in particular, by setting 



f(e,x) = fW(e,x) = Jx 2 /2e, 



(49) 



which corresponds to suppressing the noise in the inter- 
val < t < e, and afterwards taking the limit e — ► 0. 
Naturally, the free initial condition also can be written 
in the form @8J| but with f {0) (t,x) = 0. 

Quite remarkably, initial condition (|49[) is compatible 
with the structure of the solution constructed in Sec. IIVI 
and in terms of functions A(t) and B(t) corresponds to 



A(s) = 0, B{e) = -l/2e. 



(50) 
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Substitution of Eqs. flU}, © and (g7]) into Eq. |28bl) 
allows one to establish that for e <C £*, the condition 
5(e) = — l/2e corresponds to 



to — L* ~ e/3 . 



(51) 



Exactly like in the case of free initial condition (see Sec. 
W} , we have to assume that at \x\ > X(t) the dependences 
(|18a|) and (|18b|) are replaced, respectively, by Eq. (f32|) 
and V(t, x) = 0. The compatibility with condition (TT4")) is 
achieved then when the interval — X(t) < x < X(t) where 
the potential is nonvanishing shrinks to a point, the con- 
dition for which is given by Eq. (|38p . A comparison of 
Eq. with Eq. (f5"Tj) allows one to conclude that for 
initial condition (|50[) L* L/2 — e/Q and to ~ L/2 + e/6, 
which after taking the limit e — ► gives 



L/2 , t - L/2 . 



(52) 



This unambiguously defines the form of the solution for 
the case of fixed initial condition. 

In this solution, the configuration of V(t, x) is fully 
symmetric not only with respect to the change of the 
sign of x but also with respect to replacement 



t=> L-t. 



(53) 



The origin of this property is quite clear. In terms of an 
elastic string, the problem we are analyzing now is fully 
symmetric with respect to replacement (|53p . therefore 
it is quite natural that the spacial distribution of the 
potential in the optimal fluctuation also has to have this 
symmetry. 

Since we are considering the limit of zero temperature 
when the free energy of a string is reduced to its en- 
ergy, which in its turn is just the sum of the energies 
of the two halves of the string, the form of the poten- 
tial V(t, x) in the symmetric optimal fluctuation can be 
found separately for each of the two halves after imposing 
free boundary condition at t = L/2. This form can be de- 
scribed by Eqs. (2?|, ([28c]) and ([2§d| with L* = t = L/2, 
where the values of Co and Do can be obtained from Eqs. 
P0)) and (|3T)|). respectively, by replacement 



F^F/2, L^L/2. 



(54) 



The value of the action corresponding to the optimal 
fluctuation can be then found by making the same re- 
placement in Eq. (|45[) and multiplying the result by the 
factor of 2, 

S Rx {F, L) = 2S free (F/2, L/2) = ^S bee (F, L) . (55) 

Naturally, the conditions for the applicability of Eq. ([55]) 
are the same as for Eq. (|45|) (see the two last paragraphs 
of Sec. [V}. The claim that the optimal fluctuation is sym- 
metric with respect to replacement (1531) and therefore 
both halves of the string make equal contributions to its 



energy can be additionally confirmed by noting that the 
sum 



S fTee (F', L/2) + S bee {F - F', L/2) 



(56) 



is minimal when F' = F — F' = F/2. 

Like in the case of free initial condition, the form of the 
optimal fluctuation is such that the whole family of ex- 
tremal string's configurations satisfying Eq. (|42|) is char- 
acterized by the same value of energy, E = H{x(t)} = F. 
Formally, this family again can be described by Eq. (14"3")) 
where xq now should be understood not as x(0) but more 
generally as x(to). 



VII. GENERALIZATION TO OTHER 
DIMENSIONALITIES 

The same approach can be applied in the situation 
when polymer's displacement is not a scalar quantity but 
a (i-dimensional vector x. In such a case, the expressions 
for the action and for the saddle-point equations retain 
their form, where now operator V should be understood 
as vector gradient. A spherically symmetric solution of 
Eqs. (p~5|) can be then again found in the form (I L8[) with 



For arbitrary d substitution of Eqs. ifTg)) into Eqs. (TTJ 
reproduces Eqs. (|19[) in exactly the same form, whereas 
Eqs. (|2TJ)) are replaced by 



C = 2dBC, 

D = (4: + 2d)BD. 



(57a) 
(57b) 



A general solution of Eqs. (fT9| and (|57|) can be then 
written as 



A(t) = A +sign(t-t Q ) 



2(dD Q y/z 



B(t) = sign(t-to) 



D 1 



k2+d 



1/2 



C(t) 
D(t) 



where 



with 



D,/cp 2+d 



4> = 4>{t) = * 



t-tf 



La 



I+(0,d) 

2(dD y/ 2 



(58a) 

(58b) 

(58c) 
(58d) 



(59) 



(60) 



and <fr(r]) is an even function of its argument implicitly 
defined in the interval — 1 < rj < 1 by equation 



I + ($,d))=I+(0,d)\r)\ 



(61) 
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Here, I±((j),d) stands for the integral 



j±(0,d) 



dq 



l/d-l±l/2 



(l- g )V2 . 



(62) 



in accordance with which 7±(0, d) is given by the Eulcr 
beta function 



l±(0,d) 



,1 1 1 



r(|)r(i±|) 



± 



(63) 



From the form of Eqs. (|61[) and (f6"2")l . it is clear that with 
the increase of \rj\ from to 1, the function $(77) mono- 
tonically decreases from 1 to 0. It is not hard to check 
that at d = 1, Eqs. (dHD and (HI]) are reduced to Eqs. 
and (|25p. respectively. 

Exactly like in the case d = 1, for free initial condition 
one gets to = 0, L* = L, Aq = and A(L) = F/J, from 
where 



Co 



2-d F 
2 JL 



A, = 



1 



J+(0,d) 



n 2 



2L 



On the other hand, Eq. (|44j) is replaced by 



Sfree — 



4ft d 



l+d/2 



JF 



d(d + 2)(d + A) D d/2 Uq 



(64) 



(65) 



where n d = 2ir d / 2 /T(d/2) is the area of a d-dimensional 
sphere. Substitution of Eqs. into Eq. (|63|) then gives 



S ilee {F,L) = K B 



p2+d/2 



UnJ^L 1 -^ 2 



with 



K, = 



8(2-d) 1+rf / 2 (2rf) d / 2 - 1 
(d + 2)(d + 4)r(d/2) 



r(i/d + i) 



r(i/d+ 1/2) 



(66a) 



(66b) 



Naturally, at d = 1 numerical coefficient coincides 
with coefficient X in Eq. ([45]) . Like in the case of 
d = 1, the value of the action for fixed initial condition, 
x(i = 0) = 0, can be found by making in Eq. (|66a[) re- 
placement (|54l) and multiplying the result by the factor 
of 2, which gives 



S fix (F, L) = 2S frcc (F/2, L/2) = - S free (F, L) . 



The exponents entering Eq. (|66a|) and determining the 
dependence of Sfrcc on the parameters of the system have 
been earlier found in Ref. Il4| from the scaling arguments 
based on the assumption that for large L, the form of 
the optimal fluctuation involves a single relevant charac- 
teristic length scale with the dimension of x which alge- 
braically depends on the parameters of the s yste m (in- 



cluding L) and grows with the increase of L [27[. How- 



ever, the analysis of this section reveals that this length 



scale, Xq — {Cq/ Dq) 1 / 2 , tends to zero when d approaches 
2 from below, as well as the value of the action given by 
Eqs. PJ. 

This provides one more evidence that at d > 2, the 
problem with purely delta-functional correlations of a 
random potential becomes ill defined [2!| and has to be 
regularized in some way, for example, by introducing a 
finite correlation length for the random potential corre- 
lator. In such a situation, the geometrical size of the op- 
timal fluctuation is determined by this correlation length 
[l4| and its shape is not universal, that is, depends on the 
particular form of the random potential correlator. Thus, 
the range of the applicability of Eqs. (|66|) is restricted to 
< d < 2 and includes only one physical dimension, 
d= 1. 



VIII. CONCLUSION 

In the current work, we have investigated the form of 
Pl(F), the distribution function of the free energy of 
an elastic string with length L subject to the action of 
a random potential with a Gaussian distribution. This 
has been done in the framework of the continuous model 
traditionally used for the description of such systems, 
Eq. ((T|). Our attention has been focused on the far-right 
tail of Pl(F), that is on the probability of a very large 
positive fluctuation of free energy F in the regime when 
this probability is determined by the probability of the 
most optimal fluctuation of a random potential leading 
to the given value of F. 

We have constructed the exact solution of the non- 
linear saddle-point equations describing the asymptotic 
form of the optimal fluctuation in the limit of large F 
when this form becomes independent of temperature. 
This has allowed us to find not only the scaling from of 
S(F) = —Iu[Pl{F)] but also the value of the numerical 
coefficient in the asymptotic expression for S(F). 

The solution of the problem has been obtained for two 
different types of boundary conditions (corresponding to 
fixing either one or both end points of a string) and for an 
arbitrary dimension of the imbedding space 1 + d with d 
from the interval < d < 2 (d being the dimension of 
the displacement vector). Quite remarkably, in both cases 
the asymptotic expressions for S(F), Eqs. (f6"6"|) and (|6"T)) . 
are rather universal. In addition to being independent of 
temperature, they are applicable not only in the case of 
(^-correlated random potential explicitly studied in this 
work, but also (for a sufficiently large L) in the case of 
potential whose correlations are characterized by a finite 
correlation radius. Note that our results cannot be com- 
pared to those of Brunet and Derrida [3(| because these 
authors have considered a very specific regime when the 
transverse size of a system (with cylindrical geometry) 
scales in a particular way with its length L. 

Due to the existence of the equivalence between 
the directed polymer and KPZ problems, the distribu- 
tion function of the directed polymer problem in situation 
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when only one of the end points is fixed (and the other is 
free to fluctuate) describes also the fluctuations of height 
in the d-dimensional KPZ problem in the regime 
of nonstationary growth which have started from a flat 
configuration of the interface, L being the total time of 
the growth. The only difference is that the far-right tail 
of Pl (F) studied in this work in the traditional notation 
of the KPZ problem [3J] corresponds to the far-left tail 
of the height distribution function. In terms of the KPZ 
problem, the independence of the results on temperature 



is translated into their independence on viscosity. 
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